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We study the dynamics of multipartite entanglement under the influence of decoherence. A 
suitable generalization of concurrence reveals distinct scaling of the entanglement decay rate of 
GHZ and W states, for various environments. 
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The notion of entangled states is, since the early days 
of quantum mechanics, a key concept when it comes to 
distinguish between the quantum and the classical world. 
Besides this fundamental aspect, entanglement attracts 
considerable interest since it can be viewed as an indis- 
pensable ingredient for quantum information processing. 
Numerous experiments have recently been carried out in 
this area, in particular on the controlled generation of en- 
tanglement between many quantum systems [30 El [H Q 
- such as to accomplish fundamental scalability require- 
ments for quantum computation. 

An obstacle for the production and observation of mul- 
tipartite entanglement resides in its fragility under the 
unavoidable interaction with the environment. This chal- 
lenges the experimentalist, but also demands a proper 
theoretical description of multipartite entanglement in 
open systems. Despite recent progress in the understand- 
ing of decoherence processes in quantum systems, a sys- 
tematic characterization of the environment induced loss 
of many-particle quantum correlations is still lacking, and 
this is specifically due to the difficulties in quantifying 
multipartite entanglement. 

Even for bipartite mixed states, apart from the partic- 
ular case of two level systems where an exact solution is 
known [fj, the situation is far from being simple. Some 
of the widely used indicators of entanglement, such as 
the positive partial transpose and negativity [3, fail to 
detect certain entangled states, while some other entan- 
glement measures for mixed states require a high dimen- 
sional optimization procedure which only provides an up- 
per bound, unable to reliably distinguish entangled from 
separable states. Only recently a lower bound (to- 
gether with a numerically manageable upper bound) for 
the concurrence of mixed bipartite quantum states was 
derived. In the multipartite case, one usually has to deal 
with bipartite cuts, where the N individual constituents 
are partitioned into two (arbitrarily chosen) subgroups. 
The available entanglement measures for bipartite sys- 
tems thus become applicable, though different cuts or 
partitions may lead to different values of entanglement 
and, furthermore, the number of possible bipartite cuts 
increases rapidly with N . 

To improve on this situation, we will here scrutinize 
the effects of decoherence on multipartite entanglement 
measured by a suitable generalization of concurrence - 



sensitive to multipartite correlations. Moreover, since 
different experimental strategies to produce multipartite 
entanglement are subject to different sources of decoher- 
ence, we implement our approach for different types of en- 
vironment coupling, acting on different types of initially 
maximally entangled multipartite quantum states. This 
finally provides a versatile toolbox to assess the scaling 
of entanglement decay with the system size under very 
general conditions. 

To start with, we assume that each individual particle 
of the system interacts independently with the environ- 
ment and, therefore, undergoes a local decoherence pro- 
cess which is mediated by dissipation, noise or dephas- 
ing. Dissipativc effects arc described by the coupling of 
the system to a thermal bath at zero temperature, and 
can represent, e.g., the spontaneous decay of a two level 
atom induced by its interaction with the vacuum modes 
of the ambient electromagnetic field. Noisy dynamics are 
related to another limit of the thermal reservoir, when 
temperature tends to infinity whilst the coupling strength 
goes to zero. Dephasing corresponds to a situation where 
no energy is exchanged with the environment, but only 
phase information is lost. All these processes can be de- 
scribed in terms of the master equation 




where p is the reduced density operator of the system. 
The Lindblad operators Ck, describing the independent 
interaction of each particle k with the reservoir, are as- 
sumed to be of the same form for all particles and can be 
written, in the weak coupling regime and in the Marko- 
vian limit, as 
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where c, and r\ describe, respectively, the system- 
environment coupling operator and its strength. Fur- 
ther specializing to the case that each particle is a two 
level system (the generic scenario in ion trap quantum 
information processing), the operators Ci in J2J can be 
written in terms of the Pauli matrices and are given by 
c = cr_ , for the zero temperature reservoir, by ci = a - - 
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and c-i = <t+, for the infinite temperature environment, 
and by c = er + (7_ for dephasing, with I\ = T in all cases. 

It is clear that, under the action of any of the above 
environments, any initially entangled state will asymp- 
totically evolve into a separable state. We are interested 
in the rate at which multipartite entanglement decays, 
and therefore need a suitable measure of multipartite en- 
tanglement which is evaluable for general states. For this 
purpose, we use a generalization of the concurrence of 
pure bipartite states 0, C 2 (#) = ^2((^\^) 2 - Tr^), 
where the reduced density matrix g r = Tr p | (vp | is ob- 
tained as a partial trace of the bipartite state. For TV- 
partite systems, one can construct 2^ — 2 different re- 
duced density matrices, half of which in general have a 
different mixing. A generalization of C2 is therefore not 
unique, providing several inequivalent alternatives. For 
our present purpose, we will focus on the specific form 



/(2^-2)<W-^Tr^ . (3) 



where a labels all different reduced density matrices, 
i.e., there are A!/ (A — n)\n\ different terms when g a 
is obtained by tracing over n different subsystems. Cm 
vanishes exactly if \^n) is TV-separable, i.e. \^n) = 
0i=i where the are pure states of the indi- 
vidual subsystems, and Cjv adopts its maximal value 
for GHZ-states \i . . .i)j\[2. Hence, the crucial merit 
of our specific definition J3J) is that CV can account for 
real multi-partite correlations 01 ~ as opposed to other 
multipartite generalizations of concurrence 0,^3 which 
extract only bipartite correlations between single subsys- 
tems and the remainder. Also note that eq. © satisfies 
another important requirement for our present purpose - 
it allows to compare the degree of entanglement of mul- 
tipartite systems with different numbers of constituents: 
Consider an A-partite state \^>n) that factorizes into an 
A — f-partite state |^jv-i) and a one component state 
|$), \^ N ) = Itpjv-],) ® |$). In this case, the A-partite 
concurrence Cn(^n) simply reduces to the A — I-partite 
concurrence of Cn-i(^ n-i), an d thus allows for a mean- 
ingful comparison of A — I with A-particle states. 

To monitor the time evolution of entanglement, we still 
need the generalization of eq. J2J for mixed states, given 
as 



C N (p) =iniJ2 Pi C N (*i) 
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where the infimum is to be found among all sets of 
probabilities pi and pure states such that p — 

^2iPi\^ i\- I n principle, this defines an optimization 
problem that is very difficult to solve exactly. Though, 
lower bounds for Cjv(p) are available ||, and can be de- 
termined purely algebraical ly i n the regime where the 
mixing of p is moderate |l3j,[l4|. Moreover, for states of 



rank two - that we shall encounter for specific time evolu- 
tions considered further down - eq. J1J can be condensed 
into an expression equivalent to the one for bipartite two- 
level systems 0, 0] , and the required optimization can be 
performed exactly. 

With the definition (0J at hand, we can monitor the 
multipartite concurrence Cn for the solutions p(t) of 
eq. Q). We will consider two types of initial states: the 
GHZ state, |^tv)ghz = (|00 . . . 0) + |11 . . . I ))/V2, and 
the W state, \V N )w = (|00 ... 01) + |00 . . . 10) + . . . + 
1 10 . . . 00)) /VN - which are known to bear incompatible 
multipartite correlations, in the sense that they cannot 
be transformed into each other by local unitary transfor- 
mations amended by classical communication [lfij . Also 
note that these states have recently been produced 0, 
for A = 3, 4 in the lab, and that it is now within exper- 
imental reach to monitor the time dependence of their 
degree of entanglement by means of quantum state to- 
mography. 

As a first observation, illustrated in Fig. ^ we find an 
essentially perfect monoexponential decay of concurrence 
in all analyzed cases, with a decay rate 7 which depends 
on the initial condition and environment model. Only for 
an infinite temperature reservoir does C^{t) vanish after 
a finite time i sep , as illustrated in Fig. ^ consistent with 



the results obtained in |lfiill7j for a depolarizing channel. 
In contrast, for zero temperature as well as for dephas- 
ing environments does Cjv(£) only vanish for t — > 00. In 
these two cases the concurrence of the decaying W-state 
gets larger than that of the decaying GHZ-state after a 
short time. In contrast to this, the situation of the ini- 
tially prepared states, where the concurrence of the GHZ- 
state is larger than that of the W-state, is preserved for 
the infinite temperature reservoir. To allow for a con- 
sistent comparison of the decay of entanglement in arbi- 
trary environments, and, furthermore, since experimen- 
talists seek to maximize and to preserve entanglement 
by preparing pure entangled states and minimizing deco- 
herence, rather than to await its fading out, we therefore 
choose 7 rather than t scp as the figure of merit in our 
analysis of entanglement decay for variable system size. 

Consistently with the experimental scenario, we as- 
sume that the initial multipartite states are pure (or 
quasi pure, i.e. the experimental preparation of the initial 
quantum state succeeds with essentially perfect fidelity) 
0, 0- Thus, on short time scales, where the degree of 
mixing of p is expected to be small, the use of our quasi- 
pure approximation is justified. Moreover, we ex- 
plicitely verified that, even for considerably mixed states, 
this approximation provides an excellent estimate of CV 
in eq. Q with significantly reduced computational effort. 
Since we always evaluate lower bounds on concurrence, 
all decay rates shown hereafter present an upper bound 
of the actual rate at which multipartite entanglement is 
lost due to contact with the environment. 

Because of their primordial experimental interest, we 
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FIG. 1: Time dependence of the multipartite concurrence 
Cjv(t), for a N = 3-partite system initially prepared in a GHZ 
(dashed lines) or W (solid lines) state. Only for an infinite 
temperature reservoir (squares) does the state turn separable 
after a finite time i scp . In contrast, the zero temperature (cir- 
cles) as well as the dephasing environment (triangles) induce 
separability only in the limit t -> oo. In all cases, the numer- 
ical results are very well fitted by an exponential (solid and 
dashed lines). 
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focus on the multipartite entanglement of GHZ and W 
states of variable size. Figure [21 shows the scaling of their 
entanglement decay rates 7 with N, under decoherence 
induced by zero temperature, infinite temperature, and 
dephasing environments. We see that the GHZ state (top 
panel of the figure) decays into a separable state with a 
rate which increases linearly with N, except for the small- 
TV behavior of 7 for the zero temperature environment. 
Indeed, the special case of a dephasing reservoir, where 
the density matrix is always a mixture of two pure states, 
and, hence, is of rank two, can be treated analytically. 
In this particular case, the concurrence just follows the 
behavior of the two non- vanishing non-diagonal elements 
of p which uniformly decay as e~ Nrt / 2 . 

Remarkably, the situation changes quite drastically for 
the W states (bottom plot of Fig- El - I n this case, only 
the infinite temperature environment gives rise to an al- 
most linear increase of 7 with N, slightly faster than 
for the GHZ states. In contrast, for dephasing and zero 
temperature reservoirs, the decay of the concurrence is 
independent of N. Moreover, the zero temperature case 
once again allows for an analytic solution (as above, the 
rank of the state reduces to two) for all N, leading to 
Cjv(t) ~ e~ rt . Consequently, the multipartite quan- 
tum correlations of W states clearly outperform those 
of GHZ states in terms of their robustness. One might 
be tempted to attribute this to the smaller initial con- 
currence of W as compared to GHZ states. Though, the 



FIG. 2: Decay rates 7 (in units of the reservoir rate F) for 
GHZ (top) and W (bottom) states, as a function of the sys- 
tem size (i.e., the particle number) N. The different environ- 
ment models are represented by circles connected by a dot- 
ted line (zero temperature), by squares connected by a solid 
line (infinite temperature), and by triangles connected by a 
long-dashed line (dephasing), respectively. Whilst for GHZ 
states the decay rates increase roughly linearly with TV, inde- 
pendently of the specific environment, the W states exhibit 
increasing decay with system size only for the infinite tem- 
perature environment. Remarkably, the decay rate of the W 
states is size-independent for dephasing and zero temperature 
environments! 



ratio 

C N ^ GHZ ) _ J N 
C N (^w) V 'N-1 

with maximum for N — 5, approaches unity for large N. 

To summarize, we have shown that an efficient mon- 
itoring of multipartite entanglement under arbitrary 
(Markovian) environment coupling is possible, for differ- 
ent classes of quasi pure initial states typically dealt with 
in state-of-the-art experiments. Furthermore, our finding 
that GHZ states are significantly more fragile under envi- 
ronment coupling than W states might indicate a robust 
pathway to scalable quantum information processing. 
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